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Abstract. A model for axisymmetric magnetized accretion flow is proposed. The dominant mechanism of energy dissipation is 
assumed to be the magnetic diffusivity due to turbulence in the accretion flow. In analogy to the advection-dominated accretion 
flow (ADAF) solutions, a constant fraction of the resistively dissipated energy is stored in the accreting gas and the rest is 
radiated. The general self-similar solutions which describe a resistive and nonrotating flow with purely poloidal magnetic field 
is first introduced. The radial dependence of physical quantities is identical to viscous ADAF solutions. Although the main 
focus of this study is nonrotating magnetized accretion flow, for rotating flow with both poloidal and toroidal components of 
magnetic field we find the radial scaling of solutions similar to the nonrotating case. We show that the accretion and the rotation 
velocities are both below the Keplerian rate, irrespective of the amount of cooling. We show that the set of the equations is 
reduced to a second order differential equation for a nonrotating flow. The geometrical shape of the disk changes depending 
on the fraction of the resistively dissipated energy which is stored in the accreting gas. However, there is a hot low-density gas 
above the disk in almost all cases. The net accretion rate is calculated for a set of illustrative parameters. 
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1. Introduction 

The importance of the accretion processes has been recognized 
in connection with the physics of many astrophysical systems 
such as protosteller disks (e.g., Mineshige & Umemura 1997; 
Tsuribe 1999; Calvet, Hartmann & Strom 2000), disks around 
black holes or compact objects or disks at center of the active 
galactic nuclei (e.g., Lynden-Bell & Rees 1971; Kato, Fukue 
& Mineshige 1998). According to their geometrical shapes, ac- 
cretion disks are generally divided into thin and thick disks. A 
mechanism for transport of angular momentum is another key 
ingredient in theory of accretion processes and many theoreti- 
cal uncertainties remain about its nature. However, the theory 
of thin disks is well understood based on a pioneer work by 
Shakura & Sunyaev (1973). In this theory it is assumed that 
there is efficient radiation cooling in the flow. 

During recent years another type of accretion disks has 
been studied, in which the energy released through viscous pro- 
cesses in the disk may be trapped within the accreting gas (see, 
e.g., Ichimaru 1977; Narayan & Yi 1995; hereafter NY). In this 
kind of flow which is known as adection-dominated accretion 
flow (ADAF), the accreting gas has a very low density and is 
unable to cool efficiently. A full analysis of the dynamics of the 
solution was presented in an important paper by NY. Since both 
the angular momentum and the energy dissipation in the flow 
is undertaken by the turbulent viscosity similar to the standard 
a model, we can call such flows "viscous" ADAF models. 
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There are growing evidences that the accretion flows 
in many astrophysical systems involve magnetized plasma. 
Accretion disks containing magnetic fields have been studied 
by many authors (e.g., Ogilvie 1997; Hawley 2001; Casse & 
Keppens 2002). However, the standard viscous ADAF model 
completely ignores the presence of an ordered global magnetic 
field and it is assumed the turbulent magnetic field contribute a 
constant fraction of the total pressure (NY). However, we can 
expect accretion flows, in which angular momentum is trans- 
ported by a global magnetic field and the energy is released by 
Joule heating. This possibility has been addressed by Kaburaki 
(2000). He presented a set of analytical solutions for a fully 
advective accretion flow in a global magnetic field and the con- 
ductivity is assumed to be constant for simplicity. 

As for the accretion flows in presence of magnetic field, 
Schwartzman (1971) was the first to point out the importance 
of the magnetic field in an accretion process. He proposed a 
hypothesis of equipartition between the magnetic and kinetic 
energy densities and this picture as mentioned above is usu- 
ally accepted in modern picture of viscous ADAF models (e.g., 
NY). More accurate study was done by Bisnovatyi-Kogan & 
Ruzmaikin (1974), where exact nonstationary solution for a 
field amplification in the radial accretion flow was also ob- 
tained. Also, Bisnovatyi-Kogan and Lovelace (2000) suggested 
that recent papers discussing ADAF as a possible solution for 
astrophysical accretion should be treated with caution, particu- 
larly as for ignorance of the magnetic field. While they obtained 
a solution for time-averaged magnetic field in a quasispherical 
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accretion flow, an analysis of energy dissipation and equipar- 
tition between magnetic and flow energies has been presented 
(Bisnovatyi-Kogan and Lovelace 2000). In our study, we fol- 
low similar ideas in order to study magnetized accretion flow 
in which magnetic energy dissipation is stored in the flow. 

Dynamics of radiatively inefficient accretion flows has 
been the subject of many studies during recent years (e.g., 
Igumenshchev & Abramowicz 1999, 2000; Stone, Pringle & 
Begelman 1999; McKinney & Gammie 2002). In these vis- 
cous accretion models, the behaviour of the flow depends on 
the standard dimensionless viscosity parameter a. On the other 
hand, as rotation is allowed, such flows become subject to 
strong convection with completely different flow pattern and ra- 
diative properties (Igumenshchev & Abramowicz 1999; Stone, 
Pringle & Begelman 1999). Numerical MHD simulations of 
radiatively inefficient flows have been done recently by many 
authors (e.g., Hawley 2001 ; Machida, Matsumoto & Mineshige 
2001; Casse & Keppens 2002). But most of these studied have 
neglected the resistive terms in the MHD equations or the re- 
sistivity has been considered only in the induction equation 
without accounting the corresponding dissipation in the energy 
equation. 

In this paper, we want to explore how the structure of a 
steady-state thick disk depends on its resistivity by solving the 
two-dimensional MHD equations for disks which have the sim- 
plifying feature that they are self-similar in radius. We will con- 
sider the general problem of accretion flow where the field an- 
nihilation is approximated by a finite "turbulent" conductivity 
and the viscosity of the fluid is completely neglected in or- 
der to construct a fully resistive model for accretion within a 
global magnetic field. We restrict our attention to the nonrotat- 
ing accretion flow, in which the flow contains a purely poloidal 
magnetic field, and the general case will be reported in future. 
However, we will find the radial scaling of self-similar solu- 
tions for the general case where the fluid is rotating and it con- 
tains both poloidal and toroidal components of the magnetic 
field. 

This paper is organized as follows. In section 2 the general 
problem of constructing a model for resistively accretion flow 
is defined. The self-similar solutions are presented in section 
3, and the effects of the input parameters are examined. The 
implications of this analysis are summarized in section 4. 

2. Formulation of the Problem 

As stated in the Introduction, we are interested in constructing 
a model for describing magnetized accretion disks. The macro- 
scopic behavior of such flows can be studied by MHD approx- 
imation. For simplicity, the self-gravity and general relativistic 
effects have been neglected. The first assumption is true as long 
as the total mass of the disk is small compared with the product 
of the mass of central object and the angular thickness of the 
disk. Thus, the basic equations of our problem are written as 
follows: 



dt 



V x (v x B - 77J), 



dp 

| + V.(pv) = 0, 

p{% + (v.V)v) = -Vp - P W +ijxB, 

at An 



p d p 

(-) +pV.y = gdiss - gcooi 



y - 1 dt p 
V.B = 0, 



(3) 

(4) 
(5) 



(1) 



(2) 



where p, v, p, Y are the density, the velocity, the pressure and 
the gravitational potential due to a central object, respectively. 
We assume that the mass of central object is M and so, the 
gravitational potential becomes *P = -GM/r where r is the 
spherical radial coordinate. Also, J = V x B is the current 
density and r] represents the magnetic diffusivity. The term on 
the right hand side of the energy equation, gdiss, is the rate of 
heating of the gas by the dissipation and Q coo i represents the 
energy loss through radiative cooling and y is the adiapatic in- 
dex. We assume that the difference between energy dissipation 
and radiative cooling is a fraction / of energy dissipation, i.e. 
Gdiss - Gcooi = /Gdiss (e.g., NY). The parameter / measures the 
degree the flow is advection-dominated and in general it may 
depends on the details of the heating and cooling mechanisms 
and will vary with both r and 6. However, we assume a constant 
/ for simplicity. Clearly, the case / = 1 corresponds to the ex- 
treme limit of no radiative cooling and in the limit of efficient 
cooling, we have / = 0. 

We assume that the disk is stationary and axisymmetric. 
There are numerous studies, in which magnetized equilibrium 
structures have been investigated via analytical or numeri- 
cal methods (e.g., Igumenshchev & Abramowicz 1999, 2000; 
Stone, Pringle & Begelman 1999; Shalybkov & Rdiger 2000; 
Ogilvie & Livio 2001; McKinney & Gammie 2002; Rdiger & 
Shalybkov 2002). This paper will be restricted to nonrotating 
accretion flow under purely poloidal magnetic fields configura- 
tion. The whole magnetic field is maintained by the electric cur- 
rent induced in the accretion region, and the solutions are deter- 
mined only by the assumption of self-similarity. According to 
the anti-dynamo theorem (e.g., Cowling 1981), such a configu- 
ration cannot be maintained in a stationary state by an axisym- 
metric flow in the presence of the resistive dissipation. We can 
resolve this criticisms in two forms: For many problems, the 
magnetic dissipation time is very long, much longer than the 
age of the system. For such problems, one assumes the field to 
have been generated in the distant past, and restricts one's at- 
tention to its effect on current dynamics, equilibrium, stability, 
etc. 

For other problems, the dissipation time is not longer than 
the age of the system. In such cases (as in this study), one is 
making some implicit assumption about the regeneration of 
the field by (dynamo) processes outside of the computational 
regime of the problem at hand (e.g., on a 3-D time-dependent 
microscale, or beyond the boundaries of the formal computa- 
tion). Whether these assumptions are consistent with the com- 
putational details is, of course, a function of how good is the 
author's physical intuition, since such assumptions lie beyond 
the physical and mathematical equations that are actually posed 
and solved. 

There many studies, in which the time-averaged physi- 
cal quantities have been analyzed so that one can still talk 
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about stationary state and this approach has been widely 
used in literature (e.g., Lovelace, Wang & Sulkanen 1987; 
Lovelace, Romanova & Newman 1994; Shalybkov & Rdiger 
2000; Ogilvie & Livio 2001; Rdiger & Shalybkov 2002). In 
this study the physical quantities such as velocity v, magnetic 
field B etc. denote the time-averaged values and in general, 
the instantaneous value of each physical variable consists of 
a stationary part and a fluctuating component. The influence 
of the fluctuating parts is included crudely, by modifying the 
transport coefficients which would apply in the absence of 
fluctuations. Thus, we neglect the time average of fluctuating 
parts of velocity multiply magnetic field and the generation 
of magnetic field in the disk by dynamo activity is not in- 
cluded (e.g., Pudritz 1981). Since we are interested in nonro- 
tating quasispherical accretion, one may ask about the instan- 
taneous magnetic field lines. Recently, Bisnovatyi-Kogan and 
Lovelace (2000) proposed a model for magnetized quasispher- 
ical accretion, in which a sketch of the instantaneous poloidal 
magnetic field lines and the time-averaged flow velocity has 
been shown in Figure 1 of their work. As for the accretion 
flow within global poloidal magnetic field, we are assuming 
the same configuration. 

We work in spherical polar coordinate (r, 8, ip) with the ori- 
gin on the central object. We assume that the net mass accretion 



rate M ■ 



I 



2nr sin 8pv r d8 is independent of radius r. Since 



r ■ pv T is independent of r and the disk is stationary and axisym- 



metric, the continuity equation Q shows vg 
form 

l|-(pv r r 2 ) = 0. 
r z or 



(6) 



The condition of constant M is similar to what has been used 
by NY, however, is relaxed by some authors and they find so- 
lutions with outflows (e.g., Henriksen & Valls-Gaboud 1994). 
In this paper, we consider the simple case of constant net mass 
accretion rate at all radii. We can write the components of equa- 
tion of motion (equation 0) in spherical coordinates: 



dv r v z dp GM 1 
Or r Or r l An 



-py 2 cote = ~ + £;(^ B r - J t B v ), 

pv r -^ + p^ = — (J t B e - J e B r ). 
or r An 

Also, the three components of induction equation are: 

Q 

— [ r sm8(VrB g -TjJy)] =0, 
08 



-[rsin0(v r fi e -77^)] =0, 
or 

d d d 

— {rv v B T - rv r B^) + —(v^Bg) - —(rjrjg) 
or 08 Or 



+ Te (nJr) = o. 



(8) 
(9) 

(10) 

(11) 



(12) 



In equations (UOL (II It and (I12l i. the terms without 77 represent 
the advection of the magnetic field by the flow motion, while 
the terms containing 77, account for the diffusion of the field. 

From equations (II Oi and il It one can simply deduce that 
r sin 8(v r Bg - rjJ v ) = 0, or 



v r B, 



0. 



(13) 



This important equation can be understood easily in this way: 
In a stationary, axisymmetric configuration, the electric field 
can be written as a gradient of a scalar potential, and its <p com- 
ponent (i.e., E v ) which is in proportion to v r Bg - rjj^ in our no- 
tation, should vanish identically. This expression appears in the 
parentheses in equations (II 01 and ( II II . So, these equations re- 
duce into one equation, i.e. equation ( 1131 . as we showed above. 
We will find self-similar solutions so that satisfy into equation 
O- 

The energy equation (0} and the magnetic flux conservation 
give 



P( 



dp p dp 



(-)- VrirWG, 



y - 1 dr p p 2 dr 



diss? 



1 d 2 
--(r 2 B r ) + 
n dr 



1 



r sin 8 d8 



(B e sin 8) = 0, 



(14) 
(15) 



where in the above equations the components of the current 
density J are 



and takes the J r 



1 



r sin 8 d8 



(B^sinff), 



1 d 

Jg^——(rB ip ), 
r dr 



1 d dB T 
J^-l T (rBg)--^]. 
r dr d8 



To close the system of equations, we should specify the func- 
tional forms of the magnetic diffusivity 77 and the energy dissi- 
pation gdiss- In this study, we have completely neglected the ef- 
fects of viscous processes and the attention has been restricted 
only to the resistive dissipation. In other words, the accretion 
flow of our model is not infinitely conducting. For simplic- 
ity, we could assume the resistivity to be constant (see, e.g., 
Kaburaki 2000). But, we assume the magnetic diffusivity is due 
to turbulence in the accretion flow and it is reasonable to ex- 
press this parameter in analogy to the c-prescription of Shakura 
and Sunyaev (1973) for the turbulent viscosity, 



'7o 



'/() 



(16) 



where Qk = ~<jGM/r 3 is Keplerian angular velocity, and c s 
is the isothermal sound speed. NY applied similar functional 
form for kinematic coefficient of viscosity, i.e. v = a(p/pQ.^). 
Exactly in analogy to NY's prescription, we are using the above 
equation as for magnetic diffusivity 77. Note that magnetic dif- 
fusivity 77 is not constant and depends on the physical vari- 
ables of the flow and 77 in our self-similar solutions, as we 
will show, scales with radius as r 1 ^ 2 . The magnetic diffusiv- 
ity has the same units as kinematic viscosity and we assume 
that the magnitude of 77 is comparable to that of turbulent vis- 
cosity (Bisnovatyi-Kogan and Ruzmaikin, 1976; Parker , 1979 
), although we have neglected kinematic viscosity. This form 



4 



Mohsen Shadmehri: A Model for Quasi-Spherical Magnetized Accretion Flow 



of scaling for diffusivity has been widely used by many authors 
(e.g., Lovelace, Wang & Sulkanen 1987; Lovelace, Romanova 
& Newman 1994; Ogilvie & Livio 2001; Rdiger & Shalybkov 
2002). Due to arbitrary rjo it is not so important which velocity 
is used for the scaling. For example, in some studies Alfven ve- 
locity is used for the velocity in the above scaling (e.g., Fendt & 
Cemeljic 2002 ), while in other works still sound speed is ap- 
plied even in the magnetized case (e.g., Lovelace, Romanova 
& Newman 1994; Rdiger & Shalybkov 2002). Nevertheless, 
the sound speed, it seems, has an advantage; because it is the 
largest intrinsic velocity parameter in the turbulent plasma, and 
may be used universally, with or without magnetic field. In this 
study, following many previous studies, we use the above pa- 
rameterized form for the diffusivity. 

For the heating term, Qdiss, we may have two sources of 
dissipation. Magnetic reconnection may lead to energy release. 
Also, we can consider the viscous and resistive dissipations due 
to a turbulence cascade. In this study, only the resistive dissi- 
pation has been considered as the dominant process: Qdiss = 

Now, we have constructed our model and the main equa- 
tions of the model are equations 0, |7}, 0, 0, t!2i . Jl 3i . 
dl4i and J15I . In the next section, we will present self-similar 
solutions of these equations. 

3. Analysis 

3.1. Self-Similar Solutions 

To better understand the physics of our accretion flow in a 
global magnetic field, we seek self-similar solutions of the 
above equations. Of course, this method is familiar from its 
wide applications to the full set of equations of MHD in many 
research fields of astrophysics. As long as we are not interested 
in boundaries of the problem, such solutions can accurately de- 
scribe the behavior of the solutions in an intermediate region 
far from the radial boundaries. 

Before presenting the self-similar solutions, we introduce a 
convenient functional form for the magnetic field. Generally, 
we may have both poloidal and toroidal components of the 
magnetic field. We can define a positive-definite magnetic flux 
function <£(r, 0) from which we may obtain the poloidal com- 
ponents of the magnetic field: 



1 <D 
B = —V X (— r-r fy) + Bfa. 
2n rs\n0 



(17) 



We can easily show that by integrating the above equation 
over the circular area nr 2 at plane of — n/2, the result, i.e. 
<D(r, 7r/2), is the defined magnetic flux. This equation automat- 
ically satisfies the equations |15i . 

Writing the equations in a non-dimensional form, that is, 
scaling all the physical variables by their typical values, brings 
out the non-dimensional variables. We can simply show that 
a solution of the following form satisfies the equations of our 
model: 



p(r, 0) = po 



(r/r ) 3 / 2 ' 



p(r,S) = po- 



P(0) 



Wr ) 5 / 2 ' 
v r (r,0) = rn K WV(0), 
v„(r, 0) = rQ K (r)Q(0), 

Bo dm 1 



B t (r,0) = 



B g (r,0) = -Bo 
B v (r, 0) = B 



27rsin6» d6 (r/ro) 5 ' 4 ' 
30(<9) 1 



87rsin6»(r/r ) 5 / 4 ' 
b{0) 



(r/roF 4 ' 



(19) 

(20) 
(21) 

(22) 
(23) 
(24) 



where po, po, Bq and ro provide convenient units with which 
the equations can be written in non-dimensional forms and Qk 
denotes the Keplerian angular velocity, viz. Q.K(r) - ~<jGM/r 3 . 
Also, we obtain the following equations for the components of 
current density: 



J t (r, &) = — — ^-(b(0) sin 0) ' 



Je(r,0) = 



ro sin d0 

Bp b(0) 
4r (r/V ) 9 / 4 ' 



(25) 
(26) 
(27) 



j (r 0) = B ° r 3m - d ( 1 dm )] 1 

v 2nr Q 16 sin (9 d0 sin0 d6 (r/r ) 9 / 4 ' 

and considering equation (I16> magnetic diffusivity 77 becomes 



T] = T)0- 



PO 



R{0) >o 



(28) 



Substituting the above solution in the momentum, induc- 
tion and energy equations, we obtain a set of coupled ordinary 
differential equations. Note that the above solutions automati- 
cally satisfy continuity equation (|6j. Equations of motion Q, 
(|8j and (|9jl become 

R(~ - Q 2 ) = -AiP -R + ^{b 2 + 
2 2 1 67T 



30 30 d 1 



4tt 2 sin6» 16 sin6» d6 sin dO 



)]}, 



(29) 



~ dP A 2 1 dd> 

-Rn. 2 cot0= -A] — + —{—^ 

d0 4n 4n 2 sin 0d0 

30 d , 1 c/0 s , b d „ . 

( -) (hsn\0)\, 

16 sinfl d0 sin0d0 sin0d0 

Rm = - f 2 {h d ± + J^± (bsinff)l 
4^ 2 sine d0 4sin0d0 

Equations dl 2l i gives 

27r yz,-^(O0)-^-^ + 
d6 sin0d0 



(30) 
(31) 



(18) 'Jl™ 



770 Pb Snrjo d P d 

— + -rr- — [ n . „ -77;(bsm8)] = 0, 



2 R 3 d6 R sin 6 dO 



(32) 
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and equation (II 3i is cast into this form 
30 3 RV(p 



d 1 d<f> 



0. 



d6 sin 9 d9 16sin6> 4r] A i P sin 6 
and finally the energy equation fl!4i becomes 



(33) 



flf , 1 t/0 N ,l 

)] + 



2(y-l) 4tt 4tt 2 16sin0 dO smO d6 



1 



16 sin6»<i0 



(fo sin6»)] 2 



where 
Ai = 

A, = 



j?0 _ , CQs 

po(GM/r ) v'ok 



) 2 . 



B " = 4,r(^)\ 



po(GM/r ) 



VOK 



(34) 



(35) 



(36) 



where vok = (GM/ro) 1 ' 2 is the Keplerian velocity, and voa = 
#o/ yj^npo is Alfven velocity. 

We can define the ratio of the thermal to the magnetic pres- 
sures /3 as 



47T 2 Psm 2 9 



90 2 /16 + (d(p/d6) 2 + 4n 2 b 2 sin 2 9 



(37) 



where /3o = 



PO 



- 87T— . Equations d29t- (!34l > constitute 
B ( 2 /8tt A 2 i_ri— f 



a system of ordinary non-linear differential equations for the 
six similarity variables 0(0), b(0), R(6), P(0), V(ff) and Q(0). 
Indeed, the behaviour of the solution depends on boundary 
conditions which are supposed based on some physical as- 
sumptions such as symmetry with respect to the equatorial 
plane. Ogilvie (1997) studied equilibrium of magnetized non- 
accreting disks by assuming self-similarity in spherical radial 
coordinate. If we neglect resistivity, our solution reduces to the 
Ogilvie's solution, i.e. the model describes a non-accreting disk 
containing poloidal and toroidal components of magnetic field. 
However, in the case of resistive accretion flow, the scalings of 
quantities with spherical radius are still similar to those pre- 
sented by Ogilvie (1997) and some authors (e.g., Blandford & 
Payne 1982), i.e. p oc r~ 3 ^ 2 , p oc r~ 5/2 and B oc r~ 5 ^ 4 . 

The main focus of this study is nonrotating accretion flow. 
But an interesting feature of the solutions is worth emphasiz- 
ing. Using equation (I33> . we can rewrite equation i29\ as 



«!-£-«*, 'a,,**,*- »£L^ (3 8) 
2 2 167T 167r 2 ?7oP sin 2 9 



Since all the similarity variables are positive except for V, we 
see that the right hand side of the equation J38l > is positive for 
all 6 except for 6 = which is not evident from the equation. 
For nonrotating flow within purely poloidal magnetic field (i.e., 
£2 = and b — 0) which is our interest in this paper, if we sub- 
stitute from equations fl33i and J34l > into the last term of right 
hand side of the above equation, it is positively well-defined 
for all 6, even in the limit of vanishing 9. However, if we con- 
sider both the toroidal and the poloidal components of mag- 
netic field, we can still say that when 9 tends to zero, the last 



term of right hand side of equation fl38i even if tends to a very 
large value, it would be negative (because V is negative, the rest 
of the variables are positive) and so, the right hand side of the 
equation is positive (there is minus sign before the last term of 
right hand side). Thus, we should have 



V 2 , 
— +Q 2 < 1. 

2 



(39) 



This inequality implies rotation velocities at below Keplerain 
rate for the solutions of the model, irrespective of the amount 
of cooling or value of resistivity //o. Also, we can say as the 
value of angular velocity increases, the accretion velocity de- 
creases. In other words, in this model, the high accretion rate 
corresponds to low rotation rate. In this study, we shall restrict 
our attention to non-rotating accretion flow, i.e. Q. = 0. 

3.2. Nonrotating accretion Solutions 

What we have is a set of complicated differential equations 
which must be solved under appropriate boundary conditions. 
Although a full numerical solutions to the equations J29I- J34I 
would now be possible, it is more instructive to proceed by 
analyzing the model in some restrictive conditions such as 
nonrotating case. To derive this solution we set Q. — and 
b = 0, which the later condition means that we consider purely 
poloidal magnetic fields. However, to further simplify the prob- 
lem, we can not impose spherical symmetry, i.e. d/d9 = 0. 
Because if we impose spherical symmetry, we can not identify 
the actual magnetic field lines, nor can we account for nonra- 
dial forces. But we know that in resistive accretion flow the 
topology of magnetic field has a vital role. 

5/3 -y 



The same as NY we are defining parameters e 



y-l 



and £f — — . As we will see the parameter €f plays an important 

role in the behaviour of the solutions. Obviously, when we set 
Q = and b — 0, equations J3 1 1 and (I32> then immediately 
drop out and the rest of equations, i.e equations j29> . d30l and 
(I34i can be written as 



d 1 d0 30 | 



sin6» d6 sin Q dQ 16 sine 
dP A 2 



A,— 



30 1 ^0 =o 



(40) 



(41) 



(42) 



dQ 16?r 3 sin6» d6 sin 9 d9 16sin0 d6 

d6 sin 9 d9 16 sin 9 770A2 

and equation (I33i remains unchanged. Now, the set of equa- 
tions reduces to four equations J33> . (I40i . d41l and J42i which 
we want to solve. 

Clearly the expressions inside the brackets in equations 
J40> . fl41i and ( 1421 are similar and by substituting from equa- 
tion J33i . these three equations become 



9A 2 



flV0 2 



V 2 5 

R(l )=-A 1 P- , , . 

2 2 2567T 3 77oAi /> s in 2 



(43) 
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A dP 



3A 2 RV<t> dtp 
d6 ' 647T 3 77 Ai p sin 2 9 d6 

= 0. 



RV(p 2 1287r 3 77 e f A 2 



P 2 sin 2 9 3A 2 
From equation (1451 . we obtain 
1287rV f A 2 P 2 sin 2 e 



PV 



3A 2 



(44) 
(45) 

(46) 



and by substituting RV from the above equation into equations 
fi3l and (|44), we have 



1 c/P 2e f rf0 

p He + ~6~d9 



0. 



Fortunately, equation J481 is integrable and gives 
P = 



(47) 
(48) 

(49) 



where n = 2ef. Now, considering the above result, we can 
rewrite equations (I46> and 1471 as 



PV = - 



1287r 3 77 e f A 



1 <^ 2 "- 2 sin 2 0, 



3A 2 



, (1 4).(^) A , 



(50) 



(51) 



From these two equations, one can easily calculate V and R as 
functions of d> as follows 



V(0) = 
where 



1 - ^1 + 2a 2 2( "~ 2) sin 4 (9 



*"- 2 sin 2 



327r 2 /3 ?7oef 



(52) 



(53) 



3(5 + 3e f ) 

and by substituting V as function of <p into equations (15 0> or 
J5 II . we can find P. as function of (p. Thus, if we can find (/>, all 
physical variables are known. Equation J33i is the last equation 
which gives us as a function of 9. Considering equations ( I49> 
and ( 15 0> . we can rewrite equation d33l as 



d 1 d<p 
d9 sin 9 d9 
or. 



16 sin0 



+ 4-n%e t cf>"- 1 sin = 



+ 47r 2 yS o ef0" _1 



0, 



(54) 



(55) 



<^ 2 16(1 -f) 

where £, = cos 9. Now, we have a second order ordinary dif- 
ferential equation for (f>(9) which can be solved using a numer- 
ical method (in Appendix we can find an analytical solution 
for n — 1). One boundary condition can be specified on the 
equatorial plane, by assuming that field lines thread the equa- 
tor vertically: 



( d9 )e=n/2 -°- 



(56) 




Fig. 1. Self-similar solutions corresponding to y — 4/3, 770 = 
0.1 and j6o = 0.1. Top: magnetic flux function <p{9) as a function 
of polar angel 9. Middle: Radial velocity V(9). Bottom: density 
R(0). 



The second boundary condition can be specified by demanding 
the magnetic flux enclosed by the polar axis to vanish: 



0(0) = 0. 



(57) 



Thus, the complicated set of differential equations is re- 
duced to a boundary value problem of a second order differ- 
ential equation. We have obtained numerical solutions of the 
above equation for different values of / and 770. For our illus- 
trative parameters, we assume ro = lpc , Bq = 10~ 5 G, M = 10 6 
M ,/3o = 0.1 and M Q = Inpor 2 VGM/r = 1.2 x 10~ 4 M Q /year 
which gives po = 2 x 10~ 20 g/cm 3 . Figure^shows some partic- 
ular solutions corresponding to 770 = 0.1, y = 4/3 and different 
values of the advected energy /. For a fixed value of y, the so- 
lutions imply increasing cooling. This general behavior is the 
same as NY for viscous ADAF. 

The behaviour of the solutions, as the six panels in Figure 
^show, is very sensitive to the value of / or €f. The top pan- 
els show variations of magnetic flux function <p{9) as a function 
of 9. While this similarity magnetic flux is zero on the polar 
axis (this is a boundary condition), it has a maximum value 
on the equator. The maximum value of (f)(9) on the equator in- 
creases with decreasing /. This growing flux contrast between 
the equatorial region and the polar region implies that at a fixed 
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Fig. 2. Ratio of the thermal to the magnetic pressures for the 
solutions shown in Figure [0 

radius on the equator, the magnetic flux increases with decreas- 
ing y (or, equivalently, increasing cooling). 

The middle panels show similarity function of radial ve- 
locity V(ff) (see equation (I52» . We see that the maximum ac- 
cretion velocity is at equatorial region and on the polar axis 
there is no mass inflow. As we expected the accretion velocity 
is sub-Keplerian. Although for our representative set of param- 
eters the variation of V(ff) as a function of 9 is decreasing from 
equatorial region to polar region, we can recognize two distinct 
regions in the general profile of V(ff) depending on the value of 
the velocity. The bulk of accretion occurs from the equatorial 
plane at 9 = n/2 to a surface at 9 — Of, inside of which the 
velocity V(ff) is zero. Note that all of our solutions have well- 
behaved nonsigular V(ff) at all 9 and the radial velocity does 
not diverge at 9 — Of. Also, for the flows which extend from the 
equatorial plane to a surface at 8 — 6 S where 6f < 9 S < n/2, 
the accretion velocity is decreasing. With decreasing /, the re- 
gion with maximum accretion velocity becomes thin, i.e. # s de- 
creases. 

In the bottom panels of Figure ^ we plot the density sim- 
ilarity function R{9). The density contrast between equatorial 
and polar regions increases with decreasing /. The same as the 
magnetic flux function (p(9), as / decreases, the density grows 
and becomes concentrated toward the equatorial plane. As al- 
ready mentioned, the small / solutions correspond to efficient 
cooling for a fixed y. Thus, the R(9) profile clearly shows what 
we expect with increasing cooling, or equivalently ep. More in- 
terestingly, although the solutions for the large ef implies the 
bulk of the accretion occurs in equatorial region (i.e., like a 
thin disk configuration), in all cases there is a low density 
with higher temperature above the disk as the bottom panels 
of Figure [2 shows. While bulk of the accreting gas along the 
equatorial plane has approximately constant temperature, this 
tenuous region also has a constant higher temperature. As cool- 
ing increases, the temperatures of both regions decrease. Note 
that at tenuous region, there is no accretion. 

From these Figures we see for fixed values of y and 770, 
as the advected energy decreases and more energy radiates out 




-10 -5 5 10 

X 



Fig. 3. Representation of magnetic field lines in the meridional 
plane for a typical self-similar solution with 770 = 0.1, y = 4/3 
and/ = 0.84. 

of the system, magnetic flux function <p(0) increases. It means 
magnetic field lines can penetrate through the disk easier, when 
we have cooling. Thus, while thermal pressure decreases due to 
the cooling, the magnetic pressure increases and we have lower 
values for the ratio of these two which seems to be important in 
analyzing the dynamics of the flow. All solutions correspond to 
poloidal magnetic field configuration, in which the field lines 
bend only once when passing through the equatorial plane. 

Figure |2 shows the ratio of the thermal to magnetic pres- 
sures. We see that this ratio decreases by decreasing the amount 
of the advected energy /. Figure [3] shows magnetic field lines 
for a typical self-similar solution with 770 = 0.1 and / = 0.84. 
For other input parameters, we can obtain the same configura- 
tion. Ogilvie (1997) studied nonaccreting magnetized disk by 
considering a polytropic equation of state. However, the mag- 
netic field configuration in his study is similar to what has been 
obtained in this study for resistive accretion flow. If we neglect 
turbulent diffusivity and the energy equation, our solutions re- 
duce to Ogilvie's solutions. But, in our scenario, two main fac- 
tors in dynamics of the flow are the diffusion mechanism of 
magnetic field lines and the advected energy. For example, the 
net rate of accretion depends on the enclosed flux function <p{9). 
We can write 

M = -2M Q I 2 sm6R(&)V(0)dO. (58) 
Jo 

Considering equations ( 1421 and ( I49L we can rewrite the above 
equation as 

M = ^7!- 2 /3 n77oAiM \ </) 2 "~ 2 sin 3 9d9. (59) 
3 Jo 

Although this kind of magnetized accretion flow is different 
from Bondi (1952) accretion in various aspects, we can define 
Bondi accretion rate for the illustrative parameters as MBondi — 
4nG 2 M 2 p /c 3 s . Thus, M Bon di = 2(v OK /co s ) 3 M . Assuming Ai = 
1, we can write the net accretion rate as follows 

M = ^7r 2 /3o«77oM B ondi fV"- 2 sin 3 &/0. (60) 
3 Jo 
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Fig. 4. The net accretion rate M/MBondi as a function of the 
advected energy / for self-similar solutions shown in Figure^ 

For the flows which we have at hand, the mass accretion rate 
is at below Bondi accretion rate for wide range of input pa- 
rameters (e.g., ejf and 770). Figure |4] shows the net mass accre- 
tion rate M/Mfiondi as a function of the advected energy / cor- 
responding to / = 0.84,0.79,0.75,0.71,0.48,0.4,0.3,0.2 and 
y = 4/3 and 770 — 0.1. As the above equation shows the mass 
accretion rate changes in proportion to 770. For these parame- 
ters, we have fitted an approximate function as M/M Bon( ji ~ 
0.19 - 0.38/ + 0.19/ 2 . The radial structure of the resistive 
radiatively inefficient accretion flow is dramatically modified 
from Bondi profile, and the mass accretion rate is decreased 
significantly by increasing the advected energy /. We see for 
a fixed value of 770, the enclosed flux function <f>(8) increases 
as ef increases (i.e, cooling increases) and the net mass accre- 
tion rate increases, as well. Thus, in our model, the accretion 
rate is always suppresses comparing to the Bondi rate, although 
the precise amount depends on the cooling and the resistivity. 
There are many numerical simulations which show the infall of 
non-rotating or slowly-rotating magnetized gas with reduced 
M (e.g., Igumenshchev et al. 2000; Hawley, Balbus & Stone 
2001). However, our simple model for magnetized accretion 
flow shows this behavior clearly. 

4. Discussion and Conclusions 

In this paper, we have obtained axisymmetric magnetized self- 
similar advection-dominated flow solutions by considering an 
ordered magnetic field and a resistive dissipation of the energy 
of the flow. The radial dependence of the solutions has the same 
behaviour as ADAF solutions (NY), but the angular parts are 
obtained by solving a set of ordinary differential equations. In 
the model, the viscosity of the flow is completely neglected in 
order to make the contrast between the standard viscous and 
resistive ADAF solutions clear. We showed that the radial and 
the rotational velocities are well below the Keplerian veloc- 
ity, and this is independent of the resistivity and the amount 
of cooling. However, the main aim of this study was to study 
the nonrotating quasi-spherical magnetized flow directly just 



by solving the relevant MHD equations. The same problem has 
been addressed by many authors, under the assumption that the 
field will achieve equipartition of magnetic and kinetic ener- 
gies (see, e.g., Shvartsman 1971; NY; Bisnovatyi-Kogan and 
Lovelace 2000). Bisnovatyi-Kogan and Lovelace (2000) have 
discussed about equipartition between magnetic and kinetic en- 
ergies using time-averaged magnetic field solution. They ar- 
gued for such an equipartition half of the dissipated energy of 
the accretion flow results from the destruction of the magnetic 
field. In this paper we showed that the equipartition depends 
on the input parameters, more importantly the amount of ad- 
vected energy (ef) and the resistivity 770. Also, the unmagne- 
tized problem is described by pure hydrodynamics and it was 
solved by Bondi (1952) and has been widely applied. Note that 
our solutions differ in various aspects from Bondi (1952) spher- 
ical accretion. While the magnetic field diffusivity is the main 
factor in accretion in our model, the energy equation has been 
considered as well. In Bondi accretion, there is no mechanism 
for energy transport. Note that in pure Bondi spherical flow, a 
self-similar accretion is allowed only for a single value of y. 
However, magnetized rediatively inefficient accretion flow is 
possible for a range of values of y. The other interesting fea- 
ture is that the dynamics of such flow depends not only to mag- 
netic diffusivity 770, but also the fraction of advected energy / 
is important. 

For a nonrotating resistive accretion flow, the system of 
equations simplified to a second order differential equation. For 
given £f and 770, this equation determines the magnetic field 
and the other physical variables can be found using a set of 
obtained algebraic equations. Although, we have neglected in- 
gredients such as angular momentum and viscosity, the general 
properties of our solutions are similar to ADAF solutions. More 
importantly, the geometrical shape of the flow is determined by 
the amount of the advected energy (i.e., ef) for a given 770. While 
solutions with efficient cooling resemble to thin disks, the ra- 
diatively inefficient accretion describe nearly quasi-spherical 
flows. 

The bulk of accretion with nearly constant velocity occurs 
in a region which extends from the equatorial plane to a surface 
at G s , inside of which the accretion velocity gradually decreases 
to zero. Thus, the accretion velocity is constant in the region 
n/2 < 9 < 9 S . As £{ decreases, the size of this region increases 
and the accretion tends to a quasi-spherical flow. The temper- 
ature of this region is constant as well and is surrounded by a 
region with low-density and higher temperature. This hot gas 
corona is similar to what has been discussed in viscous ADAF 
solutions (NY). As the cooling increases and the flow tends to 
disk like configuration, the size of this corona increases. 

The ratio of the thermal to the magnetic pressure /3(8) de- 
pends on the input parameters €f and 770. In the limit of large e t , 
the value of (3(9) is lower than one for all 9. It means that the 
thermal pressure is not in equipartition with the magnetic pres- 
sure and more precisely, the magnetic pressure is stronger than 
thermal pressure in this regime. On the other hand, for ineffi- 
cient cooling flows, we found that the value of the ratio exceeds 
than one for a range of 9. Still there is no equipartition between 
the thermal and the magnetic pressures except for a surface at 
a specific angle which can be determined. 
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As for the rate of heating Qd\ ss , the resistive dissipation at 
small scales as a result of a turbulence cascade has been as- 
sumed as the dominant mechanism. However, there are other 
sources such as energy release through magnetic reconnection 
and more importantly viscous dissipation. Phenomenological 
considerations may lead us to define Qa ss as a fraction of the 
gravitational potential energy of the accreting gas. This simple 
definition allows one to consider all the dissipation processes 
without having detailed models for them. Preliminary results 
for this form of Qdiss shows that the scalings of quantities with 
radius still remains similar to our model, but the other proper- 
ties of the solutions should be studied. 

Since we completely neglected the viscous dissipation, the 
turbulent magnetic Prandtel number of our model is zero. This 
unsatisfactory aspect can be removed by considering both the 
viscosity and the resistivity simultaneously. Fortunately, the ra- 
dial scaling of the physical variables in the resistive case, as we 
showed, is similar to the viscous case. By dimensional analy- 
sis, one can simply show that self-similarity is possible as long 
as rj (and viscosity v) scales with radius as r 1 ^ 2 . Thus, it is still 
possible to find self-similar solutions in the general case. Given 
the importance of the value of the turbulent magnetic Prandtel 
number, we can construct a model for viscous-resistive accre- 
tion flow, in which this number is a free parameter. 
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Appendix A: 

Equation J55I can be solved analytically for n — 1 . Assuming 
16 



equation J55i reduces to 

d 2 y 3 y 

de i6i- f 

By changing the variables as 
1 

z= __ >y = z y 

this equation becomes 

d 2 Y dY 3 

z(1 -^ + (2 - 2 ^-T6 y = a 

The general solution of this hypergeometric equation is 
Y(z) = C.Fi 3 -, i|2|z) + C 2 F( 3 -, i|2| z )x 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



I 



dz 



[zF{\ ±|2|z)] 2 ' 



(A.5) 



where C\ and C2 are arbitrary constants and F is hypergeomet- 
ric function. 



